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It is shown that the parametric spectral statistics in the critical random matrix ensemble with 
multifractal eigenvector statistics are identical to the statistics of correlated ID fermions at finite 
temperatures. For weak multifractality the effective temperature of fictitious ID fermions is pro- 
portional to T e ff oc (1 — d n )/n <C 1, where d„ is the fractal dimension found from the n-th moment 
of inverse participation ratio. For large energy and parameter separations the fictitious fermions 
are described by the Luttinger liquid model which follows from the Calogero-Sutherland model. 
The low-temperature asymptotic form of the two-point equal-parameter spectral correlation func- 
tion is found for all energy separations and its relevance for the low temperature equal-time density 
correlations in the Calogero-Sutherland model is conjectured. 
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The spectral statistics in complex quantum systems 
are signatures of the underlying dynamics of the corre- 
sponding classical counterpart. The spectral statistics 
in chaotic and disordered systems in the limit of infinite 
dimensionless conductance g is described by the classi- 
cal random matrix theory of Wigner and Dyson |ij (WD 
statistics). The WD statistics possess a remarkable prop- 
erty of universality: it depends only on the symmetry 
class with respect to the time-reversal transformation 
7~. The three symmetry classes correspond to the lack 
of T-invariance (the unitary ensemble, (3 = 2); the T- 
invariant systems with T 2 = 1 (the orthogonal ensemble, 
(3 = 1), and the T- invariant systems with T 2 = — 1 (the 
symplectic ensemble, (3 = 4), respectively. The physi- 
cal ground behind this universality is the structureless 
eigenfunctions in the ergodic regime which implies the 
invariance of the eigenfunction statistics with respect to 
a unitary transformation of the basis. 

In real disordered metals the eigenfunctions are not 
basis-invariant. The basis-preference reaches its extreme 
form for the strongly impure metals where all eigenfunc- 
tions are localized in the coordinate space but delocalized 
in the momentum space. In this case the spectral statis- 
tics is Poissonian in the thermodynamic (TD) limit. 

For low-dimensional systems d = 1,2, where all states 
are localized in the TD limit, one can observe the smooth 
crossover from the WD to the Poisson spectral statistics 
as a function of the parameter £/L, where £ is the lo- 
calization radius and L is the system size. The depen- 
dence of the spectral correlation functions on the energy 
variable s = E/A (A is the mean level separation) is 
non-universal for finite L/£ but all of them tend to the 
Poisson limit as L/£ — > oo. 

In systems of higher dimensionality d > 2 the situa- 



tion is different because of the presence of the Anderson 
localization transition at a critical disorder W = W c . In 
the metal phase W < W c the dimensional conductance 
g(L) — ► oo as L — > oo and one obtains the WD spectral 
statistics in the TD limit. In the insulator state g(L) — > 
at L — > oo and the limiting statistics is Poissonian. How- 
ever, there is a fixed point W = W c in which the spectral 
statistics are nearly independent of L. Thus at the crit- 
ical point there exist universal spectral statistics which 
are neither WD nor Poissonian but rather a hybrid of 
both S. However, the universality of the critical spec- 
tral statistics (CSS) is somewhat limited, since it depends 
not only on the Dyson symmetry parameter (3 but also 
on the critical value of the dimensional conductance g* 
which in turn depends on the dimensionality d of the sys- 
tem H . Thus for each universality class there is a family 
of critical spectral statistics parametrized by the critical 
dimensionless conductance g*. 

The very existence of the subject of critical level statis- 
tics imposes a constraint on the possible values of the 
localization length exponent v = [d(3(g)/dlTig]~ 1 \ g=g * , 
where (3(g) = d\ng/dhiL is the scaling function. In- 
deed, for the spectral statistics to be meaningful the 
width of the critical energy window SE must be much 
larger than the mean level separation A oc 1/L d . The 
quantity SE is defined as the distance from the mobility 
edge E = E c at which the localization or correlation ra- 
dius £(5E) oc |<5i?|~ 1 ' is equal to the system size L. The 
number of critical eigenstates Af = SE/ A is proportional 
to L d ~ v . For vd > 1 this number tends to infinity in 
the TD limit L — > oo despite the width of the critical 
energy window shrinks to zero. This necessary condition 
for the existence of the critical statistics is secured by the 
famous Harris criterion vd > 2. 
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However, the critical exponent v enters not only in the 
necessary condition for the CSS but also in the correla- 
tion functions of the density of energy levels p(E). It 
has been shown in [|]J^] that there is a power-law tail in 
the critical two-level correlation function (TLCF) R(oj) = 
((p(E)p(E+uj))) that arises because of the finite-size cor- 
rection to the dimcnsionless conductance g(L u ) / g* — 1 oc 
(L^/L) 1 ^ = s~^, where L u oc < L is the length 

scale set by the energy difference ui = E — E' = sA be- 
tween two levels. The sign of this tail depends on whether 
the critical energy levels are on the metal (E < E c ) or 
on the insulator (E > E c ) side of the mobility edge, in 
the same way as for the 2D systems in the weak dereal- 
ization ([3 = 4) or weak localization {(3—1) regimes [ p"2| |. 
Clearly, this power-law tail does not reflect properties of 
the critical eigenstates but rather the behavior £ {\E— E c \) 
of the size of the space region where eigenstates show the 
critical space correlations. 

In order to study the relationship between the prop- 
erties of the critical eigenstates and the CSS in its pure 
form one should consider a system with a continuous line 
of critical points where /3(g) = 0. This case formally cor- 
responds to v = oo and the finite-size effects are absent. 

Another complication which makes ambiguous the def- 
inition of CSS is the fact that being independent of the 
system size, the spectral correlation functions depend on 
the boundary conditions [l3|-|l5|] and topology of a sys- 
tem. Therefore we will consider the system of the torus 
topology where CSS takes its 'canonical' form. In partic- 
ular, the TLCF decays exponentially in this case. 

As has been already mentioned the universality of the 
WD statistics is based on the ergodic, basis-invariant 
statistics of eigenfunctions which one may encounter in 
different physical situations. The characteristic feature of 
all critical quantum systems is the multifractal statistics 
of the critical eigenfunctions |]-^]. The simplest two- 
point correlations of the critical wave functions can be 
obtained from the renormalization group result [l0[| for 
I < r < £ < L: 

(\* E (0)\ 2n \* E (r)\ 2n ) =p< |* B (0)| 2 " > 2 (f/r) -, (1) 

where I is the short-distance cut-off of the order of the 
elastic scattering length, p = (£,/L) d is the probability 
for a reference point to be inside a localization region. 
The exponent a n = 2n(d n — d^n) + d^n — 2d rl + d is ex- 
pressed through the fractal dimensions d n defined by the 
L-dependence of the moments of inverse participation ra- 
tio: 



L d (\V E (r)\ 2n ) cx ir d - n>2. 



(2) 



At the critical point p = 1 and the correlation radius 
£ — > oo in Eq.(Q) must be replaced by the sample size L. 
Spectral statistics are related with eigenfunction correla- 
tions at different energies E and E' = E + u. If ui ^S> A 



and thus L u <C L one should substitute L u for L in the r- 
dependent term of Eq. (|l|) . In this way the multifractality 
exponents d„ enter the spectral ^-dependences [|| . 

For weak multifractality one can expect the fractal di- 
mensions to be a linear function of n which is controlled 
by only one parameter a: 



d n /d = 1 — an, 



1/3* « 1. 



(3) 



This relationship holds approximately for the Anderson 
transition in 2 + e dimensions and is fulfilled ex- 

actly for the critical eigenstate of the Dirac equation in 
the random vector-potential 1 1 1 . 

Thus for critical quantum systems with weak multi- 
fractality it is natural to expect that the spectral statis- 
tics depends on only one system-specific parameter - the 
critical conductance g*. 

In view of the expected universality, it is useful to find a 
simple one-parameter random matrix ensemble with the 
multifractal eigenfunction statistics which would play the 
same role for the critical systems as the classical RMT 
does for the ergodic systems. As a matter of fact there 
are few candidates ||lq| . However, here we focus only on 
one of them (ItJ , since for this ensemble the multifractal- 
ity of eigenstates has been rigorously proven [fl7 18|. 

Consider a Hermitean N x N matrix with the real 
(}3 = 1), complex ([3 = 2) or quaternionic ((3 = 4) en- 
tries Hij (i > j) which are independent Gaussian random 
numbers with zero mean and the variance: 



(l^l 2 ) 
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l/A i = j, 
1/2, i + 3 



J(i-3) (4) 



This model has been shown to be critical both for large 
p7[ and for small |Q values of B with the fractal dimen- 
sionality di at the center of the spectral band being: 



do = 



25, 



£<1 



(5) 



Thus the ID system with long-range hopping described 
by the matrix Hamiltonian Eq.([|) possesses the line of 
critical points B £ (0, oo), the fractal dimensionality d 2 
changing from 1 to with decreasing B. 

One can extend this matrix ID model by closing it into 
a ring and applying a 'flux' ip € [0, 1]. In this case 



Hij{<p) = Hij + H[] ] e 2«*> S9 «(i-i) 



(6) 



is a sum of two independent Gaussian random numbers 
with the variance of given by: 

rWi2 



H^r) = J(N-\i-j\ 



(7) 



For large values of B which correspond to weak mul- 
tifractality one can derive jl7) an effective field theory - 
the supersymmetric nonlinear sigma-model ]l9| - which 
describes the spectral and eigenfunction correlations of 
the critical random matrix ensemble Eq.(^|): 
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where Q is the supermatrix with Q| 
.9* = 4/35. 
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(9) 



The symmetry with respect to time reversal is encoded 
in the symmetry of in exactly the same way as for 
the diffusive sigma-modcl pll. The only difference is 
the long-range kernel U\i_j\ with the Fourier-transform 
Uf. = \k\. For a torus geometry k = 2irm/N, where m is 
an arbitrary integer. 

One can explicitly resolve the constraint Q 2 = 1 by 
switching to the integration over the 'angles' W. Then 
the Gaussian fluctuations of 'angles' recover the spectrum 
of 'quasi-diffusion' modes: 



5 *H, m = 0,±l,±2,.. 



(10) 



The problem of spectral statistics can be generalized to 
include the dependence of spectrum on the flux tp in- 
troduced by Eq.(g). One can define the parametric 
two-level correlation function R(s,<p) = ((p(E,0)p(E + 
sA,<p))) which can be treated in the framework of 
the same nonlinear sigma-model but with the phase- 
dependent e m : 



i(f) = 9* \m-cp\. 



(11) 



Following the work by Andreev and Altshuler Ref. |21 
we introduce the spectral determinant: 



D-\s,<p) 



n 



,(*>) 



4,(0) 
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Then it can be shown in the same way as in Ref. 21 that 
the parametric TLCF for s^> 1 and g*tp ^> 1 can be ex- 
pressed in terms of the spectral determinant as follows: 
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ds 2 
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d 2 G(s,tp) 


2tt 2 


ds 2 
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d 2 G(s,tp) 


8tt 2 


ds 2 



■ cos(27rs) e G{s ^ (13) 
2cos(27rs)e 2G ( s ^ (14) 



— = cos(27rs) e 
V8 



G(s, V )/2 



+ 



R u (s,v) = - 

R°(s,<p) = - 

R s (s^) = - 

+ I cos(47rs)e 2G(s ' v) , (15) 
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where G(s, ip) — G(s, tp + 1) is a periodic in ip function: 



2ir 2 (s 2 +e 2 ( l p))- 



(16) 



Eqs.(n3|jTq) coincide with the corresponding formulae in 
Ref. |2l| ] for the unitary, orthogonal and symplectic en- 
sembles after some misprints are corrected as in Ref. |22] 



and s — ► 2s for the symplectic ensemble to take account 
of the Kramers degeneracy. The only difference is in the 
form of the spectral determinant Eq. (|lj) due to the spe- 
cific spectrum e m (tp) of the quasi-diffusion modes. 

Using the functional representation Eq. (0) one can also 
find the leading term in the deviation from the WD statis- 
tics at s <C g* and tp — using the results of Ref. |^3],|l7| : 

where Rwd(s) is the Wigner-Dyson TLCF. 

A remarkable property of the function G(s, ip) for the 
critical sigma-model Eq.(^) on a torus is that it can be 
decomposed into the sum G(s, tp) = F(z) + F(z) of ana- 
lytic functions F(z) and F(z) where r = g*tp, z = r + is, 
z = t — is and 



F{z) = - ln[V2g* sm(nz/g*)}. 



(18) 



Eq.(|18[) results from a straightforward evaluation [|16[ of 
the product in Eq.(|l2|). 

On the other hand, it can be easily verified that G(s, tp) 
given by Eq. (|l8|) is proportional to the Green's function 
of the free-boson field &(z) on a torus in the (1 + 1) z- 
space: < < g*, — oo < 3z < +oo: 

(*(«, ,g» $(0, 0)) s - ($(0, 0) $(0, 0)) S = K G(s, ip), 

(19) 

where (...)s denotes the functional average with the free- 
boson action: 



i 



8irK 



g~ r+oo 

dr I ds [(<9 S $) 2 + (<9 T $) 2 ] . (20) 

oo 



Now we are in a position to make a crucial step and 
suggest that for the critical RMT described by Eq.([|), 
the Andreev- Altshuler equations ( ^3[jl5| ) are nothing but 
density-density correlations in the Luttinger liquid of fic- 
titious ID fermions at a finite temperature T = 1/g*: 

R(s,ip) = n~ 2 (n(s,T)n(0,0)) s - 1, r = g V (21) 

Indeed, the density operator n(s, t) (s- is space and 
r G (0, 1/T) is imaginary time coordinate) for ID in- 
teracting fermions with the Fermi- momentum kp = it 
can be expressed through the free boson field &(s,t) as 
follows |3 : 

n (s, T ) ~ 77~" d s®(s, t) + A K cos[27rs + $(s, r)] + (22) 
Ztt 

+ B K cos[47rs + 2$(s, t)]. 

The constants Ak and Bk are independent of 'temper- 
ature' 1/g* but depend on the interaction constant K. 
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They can be uniquely determined from the WD limit 
9* -> oo. 

Using Eqs.( ^l| , p2[ ) and the well known result for the 
Gaussian average of the exponent: 



ip<S>(s,T) -ip$(0,0)\ _ Kp 2 [G(s,t)-G(0,0)] 



one can verify that for the choice: 



K=-, (3 = 1,2,4 

P 



(23) 



(24) 



the Andreev-Altshuler formulae Eqs.(|l3[]l5|) are repro- 
duced exactly for the orthogonal, unitary and symplectic 
ensembles, respectively. Now we remind on a known re- 
sult that the parametric spectral statistics in the WD 
limit g* = oo is equivalent to the Tomonaga-Luttinger 
liquid at zero temperature. It follows directly from 
Ref. pf| and the equivalence of the Calogero-Sutherland 
model and the Tomonaga-Luttinger liquid for large dis- 
tances \z\ 3> 1. The critical random matrix ensemble 
Eq.(^,UJ^) and the critical ID sigma-model Eq.(^) turns 
out to be the simplest generalization of the WD theory 
that retains the Tomonaga- Luttinger liquid analogy ex- 
tended for finite 'temperatures' T = l/g* which are re- 
lated with the spectrum of fractal dimensions Eq. (|j) . 

The Wigner-Dyson two-level statistics for all three 
symmetry classes can be expressed through the single 
kernel K(s) = sin(7rs)/(7rs) in the following way [EJ: 



R°(s) = -K 2 (s) - 
R°(s) = -K\s) 



R u (s) 
dK(s) 
ds 

dK(s 



ds 



K(x) dx, 
K (x) dx. 



(25) 
(26) 

(27) 



It turns out that such a representation is also valid for 
the critical TLCF at g* 3> 1 if the kernel is replaced by: 



K(s) 



T sin(7ras) 
a sinh(7rTs) 



T 



1 
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(28) 



where a = 1 for the orthogonal and the unitary ensemble 
and a = 2 for the symplectic ensemble. The form of the 
kernel Eq.(28) can be guessed from the well known den- 
sity correlation function for the case of free fermions in 
one dimension at a finite temperature that corresponds 
to the unitary ensemble |l6|,^5|,|7) . In order to prove this 
statement we note that for s ^ 1 and g* ^> 1 Eqs.(p5-27) 
with the kernel Eq.(p8|) give the same leading terms as 
Eqs.(p[|l|) with G(s,0) given by Eq.(|l|). For s < g* 
Eqs. (25-27) give corrections to the WD statistics that 
coincide with Eq. (|l7|) . Thus the representation Eqs.(p5f- 
p7| ) with the kernel Eq.(j28|) are the correct asymptotic 
expressions for both s ^> 1 and s <C g* ■ At g* ^> 1 
these regions have a parametrically large overlap so that 
Eqs.(p5|]27]) are valid for all values of s. Given that for 



all s at T = and for large energy separations s 1 
at T <C 1 expressions Eqs. (^5- 27) correspond to equal- 
time correlations in the Calogero-Sutherland model, one 
can expect Eqs.(|25|-p7|) with the kernel Eq.(p8|) to de- 
scribe the low-temperature equal-time correlations in the 
Calogero-Sutherland model Q for for all s and the in- 
teraction constant K = 2/(3. 

One can use the kernel Eq.(28) to compute the spac- 
ing distribution function P(s, g*) which can be compared 
with the corresponding distribution function P c (s) ob- 
tained by the numerical diagonalization of the three- 
dimensional Anderson model at the critical point. Such 
a comparison is done in Ref. |28j for (3 = 1,2,4. It 
turns out that identifying the parameter g* from the fit- 
ting P(s,g*) ~ e~ K ( 9 ) s with the far exponential tail of 
P c (s) ~ e~ KS one reproduces the entire distribution func- 
tion P c (s) extremely well. 

A. M. T. acknowledges a kind hospitality of Abdus 
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